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ABSTRACT 

Recent measurements showed that the period derivative of the 'high-B' radio pulsar 
PSR J1734— 3333 is increasing with time. For neutron stars evolving with fallback 
disks, this rotational behavior is expected in certain phases of the long-term evolution. 
Using the same model as employed earlier to explain the evolution of anomalous X-ray 
pulsars and soft gamma-ray repeaters, we show that the period, the first and second 
period derivatives and the X-ray luminosity of this source can simultaneously acquire 
the observed values for a neutron star evolving with a fallback disk. We find that 
the required strength of the dipole field that can produce the source properties is in 
the range of 10 12 — 10 13 G on the pole of the neutron star. When the model source 
reaches the current state properties of PSR J1734— 3333, accretion onto the star has 
not started yet, allowing the source to operate as a regular radio pulsar. Our results 
imply that PSR J1734-3333 is at an age of - 3 x 10 4 - 2 x 10 5 years. Such sources 
will have properties like the X-ray dim isolated neutron stars or transient AXPs at a 
later epoch of weak accretion from the diminished fallback disk. 
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1 INTRODUCTION 

The discovery of several new classes of isolated neutron 
stars, namely the anomalous X-ray pulsars (AXPs) and soft 
gamma-ray-burst repeaters (SGRs), the X-ray dim isolated 
neutron stars (XDINs), the compact central objects in cer- 
■ tain supernova remnants (CCOs), and the rotating radio 
transients (RRATs), has brought into focus the question 
of po ssible evolutionary links (see, e.g.. Kaspi 20101 : iPopovl 
l200ct ). The existence of radio pulsars with large inferred 
dipole magnetic moments, close to, and in fact partly over- 
lapping with, the range of magnetar fields inferred for AXPs 
and SGRs (> 10 14 G), and the observation of radio pulses 
from some AXPs and SGRs, further highlight the possibility 
of links and raise questions about the similarities and differ- 
ences among these sources. The locations and evolutionary 
tracks of pulsars in the P — P diagram hold the keys to de- 
ciphering the links. The recent measurement of the peculiar 
braking index n = 0.9 ± 0.2 of PSR J1734-3333, the lowest 
among the measured braking indices of young pulsars, with 
peri od P = 1.17 s and pe riod derivative P = 2.28 xlO -12 s 
s _1 l|Espinoza et alj|201lf ) is an exciting new clue. 

Analysing the P — P diagram of all isolated pulsars 
in terms of evolution by rotationally powered dipole radia- 
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tion into the vacuum, starting with the initial rotation rate 
and the magnetic dipole moment at birth, and assuming 
that the dipole moment remains constant, fails to explain 
the distribution of all young pulsars on the P — P diagram 
or to shed light on the possible connections of the distinct 
classes. Evolution of the magnetic-dipole moment and its an- 
gle with the rotation axis is one direction of extending the 
picture, as suggested i n the work reporting t he braking index 
of PSR J1734-3333 (|Espinoza et alj|201lh . Magnetic-field 
evolution and amplification to lead to magnetar strength 
10 14 — 10 15 G surface fields h ave been invoked to explai n 
the bursts of SGRs and AXPs |Duncan fc Thompson|[l992f ). 
This does not require that the magnetar strength fields are 
necessarily in the dipole component that controls the spin- 
down of the pulsar and thereby its evolution on the P — P 
diagram. Indeed, the recently pu blished P upper limit of 
SGR 0418+5729 |Rea et alJl201Ch has shown that the sur- 
face dipole field of this pulsar is at most 7.5 x 10 12 G at 
the equator (1.5 xlO 13 G at the poles). This has been inter- 
preted in ter ms of the decay of th e magnetar dipole moment 
of this SGR |Turolla et alj|201ll ). All measured braking in- 
dices l |BeckeJ l2009) deviate from the n = 3 value, charac- 
teristic for dipole radiation in vacuum. Braking indices n < 
3 imply a growing magnetic dipole moment (or growth of 
the dipole moment component perpendicular to the rota- 
tion axis). Thus, to understand the behavior of pulsars and 
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magnetars in different parts of the P — P diagram requires 
growth, or decay, of the magnetic dipole moment perpendic- 
ular to the rotation axis, at rates depending on the sources. 

An alternative avenue for a more general picture of iso- 
lated pulsar evolution is to allow for the possibility of in- 
teraction with matter around the star, so that the emission 
is not dipole radiation in vacuum. An effective possibility is 
that some matter left over from the supernova explosion is 
actually bound to the neutron star, in the form of a 'fallback' 
accretion disk, as it necessarily carries angular momentu m, 
as proposed bv lChatteriee. Hernquist fc Naravanl (|200Cj ) to 
explain the properties of AXPs. Alparl ( 2001 ) suggested that 
the presence or absence, and the initial mass, of a fallback 
disk could be the third initial condition, complementing the 
initial rotation rate and dipole magnetic moment, to deter- 
mine the subsequent evolution of different classes of neu- 
tron stars. A first simple application of this idea to pulsar 
braking indi ces and motion across th e P — P diagram was 
presented bv lAlpar. Ankav &: Yazganl ((2001). Discoveries of 
radio pulsars with long periods and large period derivatives 
suggested that these so urces could have evolutionary links 
with AX Ps and SGRs dKaspil |2010| ; lEspinoza etafl l201ll ; 
see, e.g., lMereghettill20oj for a recent review of AXPs and 
SGRs). 

In the present paper we apply the fallback disk model 
in detail to the evolution of PSR J1734— 3333 and show 
that the model can explain all the properties of this source, 
including its braking index. In Section 2 we outline the 
model and examine evolutionary tracks for a neutron star 
with a fallback disk to search for scenarios leading to the 
present properties (P,P,P, and L K ) of PSR J1734-3333. 
We trace all possible initial conditions, namely the initial 
period, dipole field, and disk mass, that can produce the 
source properties. We discuss the results of our model cal- 
culations in Section 3, and summarise our conclusions in 
Section 4. 



2 THE MODEL 



In the fallb ack disk model (FDM) (|Chatteriee et al.ll2000l ; 
lAlparl 1200 ll ) the period and luminosity evolution are de- 
termined by the interaction between the neutron star and 
a fallback disk around it. The mid-infrared detections 
from the AXPs 4U 0142+ 61 dWang. Chakrabartv fc Kaplan! 
l200d ) and IE 2259+586 ^Kaplan et al.ll2009T) are consistent 
with the presence of fallback disks aro und these sources 
Ertan et all 120071 '). It was shown by lErtan fc Calis kan 



2006) that the observed near-IR luminosities and the up- 



per limits of AXPs/SGRs are compatible with the expec- 
tations of FDM. The values of the dipole-field strength at 
the pole of the neutron star Bp, indicated by FDM fits to 
optical/IR data of 4U 0142+61 (|Ertan et alj|2007l ) and by 
the results of the work explaining the long-term P, P, and 
X-ray luminosity, L x , evolu tion of AXPs and SGRs, are 
less t han 10 13 G in all cases (|Ertan et al.ll2009l ; lAlpar et ail 
l201ll ). It was proposed in these papers that the magnetar 
strength fields needed to power the bursts must be resid- 
ing in quadrupole and higher multipole components. The 
higher multipole fields fall off with distance from the star 
more rapidly than the dipole field does, leaving the dipole 



field to determine the interaction with the disk and the re- 
sulting torques. 

The interactive evolution of the neutron star and the 
fallback disk can have epochs with accretion as well as radio 
pulsar epochs. The neutron star enters the accretion regime 
and experiences an efficient torque if and when the inner 
edge of the fallback disk penetrates into the light cylinder. 
The neutron star can then spin down to long periods of 
several seconds on timescales from ~ 10 3 to ~ 10 s yr, de- 
pending on the disk torques, the dipole-field strength Bo 
and the disk mass Md- The first (ejector) phase of evolution 
without accretion could last from several years to more than 
10 yr depending on Bo, Ma, and initial period Pq. During 
this phase, the neutron star is a radio pulsar. In the present 
work, we show that PSR J1734— 3333 is likely to be in the 
radio pulsar phase without accretion at present and that the 
accretion epoch could start at a future time. In the accre- 
tion phase, at a time depending on the initial parameters, 
the inner disk will reach the light cylinder and the accretion 
will stop. After the accretion phase, the disk could remain 
attached to the light cylinder and the disk torque could still 
remain active while its efficiency gradually decreases to the 
level of the dipole radiation torque. Unlike a steady-state 
disk in a binary, where the accreting stage can be sustained 
on the evolutionary timescales of the binary and the com- 
panion, the fallback disk around an isolated neutron star 
will diminish. From the beginning, the outermost parts of 
the disk are always at low temperatures. Eventually, tem- 
perature even in the inner disk becomes too low to sustain 
viscosity. The disk then becomes passive, mass inflow and 
disk torques terminate. The decreasing luminosity and disk 
torque together l e ad to the observe d period clustering (see 
lErtan et al"1l2009l ; lAlpar et al.ll201ll . for details). 

For PSR J1734— 3333, we investigate the evolution 
mainly in the initial radio pulsar phase. We address the pe- 
culiar braking index of PSR J1734— 3333, thereby exploring 
the effect of a possible fallback disk on the evolution of iso- 
lated radio pulsars across the P — P diagram. The recent 
measurement of t he second derivative of the period, P = 
5.3 x 10 -24 s s -2 l|Espinoza et al.|[201 ll ). provides an oppor- 
tunity to test FDM evolutionary scenarios more stringently 
than before, checking for the first time for simultaneous 
agreement of P with P, P, and X-ray luminosity, L x . The X- 
ray luminosity of the source is 7.3 x 10 31 -6.6x 10 32 erg s , 
taking int o account the 25% e rror margins for the distance 
estimate jOlausen et al.l l2012). It is interesting that the P 
of PSR J1734— 3333 is positive. This means that the pulsar 
is evolving towards the upper right corner, the AXP/SGR 
region, of the P — P diagram. For an isolated neutron star 
evolving by magnetic-dipole radiation in vacuum, this would 
require a dipole field growing in time. Note that the toroidal 
and dipole fields of AXPs/SGRs, starting from the early 
phase of evolution, should decrea se rather rapidly wi th time 
in the magnetar model (see, e.g.. iTurolla et al"1l201lf ). 

The model we employ for PSR J1734— 3333 is the same 
as the one we used to investigate the long term evolution of 
AXPs/SGRs in our earlier work. The details of the mode l 
are described in lErtan et"aH l|2009h and lAlpar et all (|201lf ). 
We start with an initially extended disk with an inner radius 
equal to the Alfven radius 
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where M[ n is the rate of mass flow arriving at the inner 
disk, G is the gravitational constant, M and fi are the mass 
and the magnetic-dipole moment of the neutron star. When 
the inner disk radius n n , calculated by Equatio n 1, exceeds 
the l ight cylinder radius tlc, we set r- m = tlc l|Alpar et al.l 
l201ll ). This assumes that the inner disk remains linked on 
the closed field lines when it cannot enter the light cylinder. 
In this phase, accretion is not possible, and pulsed radio 
emission is allowed. In the phase of spin-down with accre- 
tion, when the inner disk is inside the light cylinder and 
greater than the co-rotation radius (r co < ri n < tlc) , a 
fraction of the matter could be propelled from the system 
while the remaining fraction accretes onto the neutron star. 

In the accretion phase, we calculate the disk torque act- 
ing on the neutron star using 

N = ~ M iB {GMr in ) 1/2 (1 - ul) = J fi, (2) 

ijErtan fe Erkutll2008h . where I is the moment of inertia of 
the neutron star. The fastness parameter is defined as lu, = 
f2»/fiK(rin), where f2K(n n ) is the angular frequency of the 
disk at ri n = va and Q,* is the angular frequency of the 
neutron star. Using Equations 1 and 2, it is found that P oc 
B 2 , independent of M and P when the system is not close to 
rotational equilibrium. This indicates that P is constant and 
P = in this phase. When a high-luminosity AXP/SGR is 
approaching (or receding from) rotational equilibrium P is 
negative (positive). In the early radio phase, the magnetic 
dipole radiation torque could dominate the disk torque for 
fast born pulsars (Po ^ 50 ms) for Bo ~ 10 12 G. We calculate 
the total torque including also the magnetic dipole torque 
N dip = ~2^ 2 ^/3c 3 . 

In our model calculations, we find that over the spin 
history of PSR J1734— 3333, the dipole radiation torque re- 
mains 2-3 orders of magnitude weaker than the disk torque 
for an initial period Po ~ 300 ms. For a given field strength, 
the ratio of the torques depends on the chosen disk torque 
model and the initial per iod Pp. For instance , our results 
agree with the results of lYan, Perna fc Sorial <|2012l ), who 
found that for slow-born pulsars (with Po ~ 300 ms) the disk 
torque dominates for the first ~ 10 5 years. In their model, 
for fast-born pulsars (Po ~ 5 ms), the dipole torque is the 
dominant mechanism for the first ~ 10 J years. The disk 
torque we employ in ou r model is more efficient than that 
used in lYan et al.l ((2012) , and depending on the disk mass, 
it could start to dominate the dipole torque in an earlier 
phase of the evolution for fast-born pulsars (which probably 
represent a small set among newborn neutron stars, lying in 
the tail of a Gaussian distribution with mean 300 ms and 
standard deviation 150 ms, acc ording to the simulations of 
iFaucher-Giguere fc Kaspi 2003). 

The mass flow rate at the inner disk is calculated at each 
evolutionary step by solving the diffusion equat ion with an 
initia lly thin disk surface density profile (e.g., lErtan et al.l 
2009). Initially, for numerical reasons, we set the outer disk 
radius at r out = 5 x 10 14 cm. After the first time step, 
Tout has a dynamical evolution such that the temperature 
at r out remains equal to T p , and beyond this radius the disk 
is assumed to be in a viscously inactive phase. Because of 
decreasing irradiation flux, temperatures and r out also de- 
crease gradually with time. The initial disk mass is found 
by integrating the initial surface density profile. Since the 



position of the initial outer disk radius is not well known, 
our Md values may not reflect the actual full disk mass. In 
this phase, we substitute r- m = tlc to calculate the disk 
torque (Equation 2). When the dipole torque is negligible, 
we find P oc Mi n P 7 ^ 2 in this phase. We perform numeri- 
cal calculations to follow the evolution of Mj n t ogether with 
corresponding P, P, an d P at each time step (|Ertan et al.l 
l2009l ; lAlpar et al]|201ll ). We repeat the calculations until we 
identify the initial conditions that can produce the observed 
P, P, P and also, in the present case, L x of PSR J1734— 3333 
simultaneously, at an age when the disk does not penetrate 
the light cylinder, allowing for radio pulsar activity. 

In the radio pulsar phase, for a neutron star with the 
properties of PSR J1734— 3333 the source of the X-ray lu- 
minosity is very likely to be the cooling luminosity. For the 
fast born pulsars, the dipole radiation luminosity could re- 
main well above the cooling luminosity until the periods 
increase to a few 100 ms. We include the cooling and dipole 
luminosities i n addition to dissipation due to magnetic and 
disk torques (|Alparl 120071 ) in the calculation of the total 
luminosity in the radio phase. In the long-term evolution 
of AXPs/SGRs and XDINs, when the sources reach long 
periods of a few seconds, the presence of cooling luminos- 
ity extends the life time of the disk by providing irradi- 
ation even in the absence of accretion luminosity. During 
the initial radio phase, it does not have a significant ef- 
fect on the evolution, but, through the effect of irradiation 
on the disk evolution and disk torques, it affects the time 
at which the model source acquires the observed rotational 
properties (P, P, P) simultaneously. In our model, we em- 
ploy the intrinsic cooling lumin osity evolution calculated by 
iPaee. Geppert fc Weberl l|2006t ) with the assumption of an 
isothermal neutron star with radius 12 km and mass 1.4 Mq. 
The results are not sensitive to the choice among standard 
cooling scenarios. 



3 RESULTS AND DISCUSSION 

The model curves that can account for the properties of PSR 
J1734— 3333 are seen in Figures 1-4. Performing many sim- 
ulations tracing the initial conditions, we obtain reasonable 
results for the range of disk masses M c \ — 9.0 x 10~ 8 — 1.5 x 
1O" 6 M and for dipole fields with So = 9 x 10 11 -1 x 10 13 G. 
The model is in very good agreement with P, P, P, and L x of 
the source with Bo — 2x 10 12 G (solid curve in Figure 1). For 
the solid curve in Figure 1, the disk mass Md = 3 x 1O _7 M0 
and the initial period Po = 300 ms. In Figure 1, we also 
present illustrative model curves (dashed and dot-dashed) 
that cannot represent the evolution of PSR J1734— 3333. 
For the viable model in Figure 1, and all three viable models 
shown in Figure 2, the source is powered by intrinsic cool- 
ing when the observed properties of PSR J1734— 3333 are 
attained. For this radio pulsar, considering all reasonable 
model curves, the simultaneous agreement with all these 
model parameters leads to a model estimate of the present 
age ~ 3 x 10 4 — 2 x 10 5 years. At an epoch later than the 
present age of PSR J1734— 3333, the disk eventually pene- 
trates the light cylinder. Accretion then starts, and the star 
enters a constant P phase. This is also observed as a small 
abrupt rise in the luminosity curve. There is no substantial 
increase in the luminosity when accretion starts, because 
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Figure 1. Evolution of the luminosity, period, first and second 
period derivatives of the model sources. Horizontal lines show the 
properties of PSR J1734— 3333, with the observational uncertain- 
ties in L x and P. The vertical lines are to show the time period 
over which the solid (red) model curve traces the uncertainty 
range of P. Values of Bo are given in the second panel. For these 
calculations, we have taken M d = 3 X 10 — 7 Mq and Pq = 300 ms. 
In the accretion phase, sources enter the constant P phase and P 
becomes (see text for details). 



Figure 2. Evolution of the luminosity, period, first and second 
period derivatives of model sources. Horizontal dotted lines repre- 
sent the properties of PSR J1734— 3333 with the range of uncer- 
tainties in L x and P. These illustrative model curves are obtained 
for Bo = 2 X 10 12 G. The values of initial period and disk mass 
are given in the second panel. The cooling luminosity is shown 
with the dot-dot-dashed (black) curve. It is seen that the source 
properties could be well reproduced with different initial periods. 
Between the vertical lines given with the same color, the model 
sources trace the uncertainty range of P (see text for details). 



the mass-flow rate has already decreased to low levels by 
the time the inner disk penetrates the light cylinder. Rea- 
sonable model curves indicate that accretion will start at a 
future time of order ~ 10 4 yr from the present. For other 
model source histories, accretion can start at earlier times 
and the accretion luminosity can be orders of magnitude 
greater than the cooling luminosity. Evolutionary curves for 
such models do not simultaneously produce all the observed 
properties of PSR J1734-3333. 

Figure 1 shows that P reaches values of order ~ 2 x 
10 -22 s s -2 at the end of the radio pulsar epoch. Before the 
accretion phase, from t ~ 10 3 years to the present age, the 
braking index for model sources varies from ~ 5 to ~ — 1. 
In the accretion phase, P — 0, P becomes constant and the 



braking index remains ~ 2 when the source is not very close 
to rotational equilibrium. 

The mass accretion from the disk starts at a time that 
depends on the in itial period Pq, as well as on Bo and M<j 
l|Ertan et al.l l2009h . A large range of Po values are allowed 
for producing the source properties. In Figure 2, for a given 
dipole field (Bo = 2 x 10 12 G), we present three illustra- 
tive model curves with initial periods of 50, 100 and 300 
ms. To obtain acceptable models, smaller Po values require 
greater Md- All model curves given in Figure 2 produce 
the X-ray luminosity and the rotational properties of PSR 
J1734— 3333 at different epochs. In the early radio-pulsar 
phase, the source of the luminosity is either the intrinsic 
cooling or magnetic dipole radiation of the neutron star, de- 
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Figure 3. Evolution of the luminosity, period, first and sec- 
ond period derivatives, for models that do not work for PSR 
J1734— 3333. Horizontal dotted lines show the properties of PSR 
J1734— 3333, with the uncertainties in L x and P. These model 
curves arc obtained with B = 2 X 10 12 G and P = 300 ms. The 
solid lines are the same as the solid lines in Figures 1 and 2, for the 
model that works. We also present two illustrative model curves 
for smaller and greater M<j that cannot represent the evolution 
of PSR J1734-3333 (see text for details). 

pending on Bo and p). Reasonable model curves indicate 
that current luminosity of PSR J1734— 3333 is produced by 
the intrinsic cooling luminosity while the source is slowing 
down by the disk torques. We find that the rates o f dissi- 
patio n inside the star due to dipole and disk torques (|Alpar| 
l2007f ) do not contribute significantly to the total luminos- 
ity in the radio-pulsar phase of this source at present. A 
part of the dipole luminosity is emitted in the X-ray band. 
For model sources that have low Po and/or high Po values, 
the dipole radiation could be the dominant luminosity in 
the early phases of evolution. For instance, for initially fast- 
rotating model sources with Po = 2 x 10 12 G, the luminosity 
of magnetic dipole radiation dominates the intrinsic cooling 
luminosity until the period increases to beyond ~ 200 ms 
(see Figure 2). For Po <; 10 13 G, the dipole luminosity re- 



Figure 4. Evolution of the luminosity, period, first and second 
period derivatives, for model sources with the lowest and highest 
allowed Bq values. The magnetic field is given in units of 10 11 G, 
the initial period is in seconds and the disk mass is given in units 
of 10 — 8 Mq. The cooling luminosity is shown with the dot-dot- 
dashed (black) curve. 

mains above the upper limit of the observed luminosity when 
the rotational properties of PSR J1734— 3333 are acquired. 
We note that the bl ackbody nature of the X-ray spectrum 
l|OIausen et al.|[2012f ) also indicates that the source of the ob- 
served luminosity is likely to be the intrinsic cooling rather 
than magnetic dipole radiation. For the model sources with 
Bo < 9 x 10 11 G, accretion starts before the current rota- 
tional properties of PSR J1734— 3333 are reached. To sum 
up, the luminosity, P, P and P values of PSR J1734— 3333 
can be acquired by the model sources only for Po values in 
~ 9 x 10 11 - 1 x 10 13 G range. The two model curves with 
the minimum and maximum allowed values of Po and Aid 
are given in Figure 4. 

The important parameters that affect the long-term 
evolution of the model sources are the initial period Po, the 
minimum temperature of the active disk T p , which is degen- 
erate with the irradiation efficiency C, ( but C is restricted 
by optical and IR observations of AXPs; lErtan fc Ca hskan 
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Figure 5. The evolution of the three model sources of Figure 
3 on the P — P diagram. The values of initial disk mass Md 
are shown in the Figure. All sources start with Po = 300 ms. The 
model source with the lowest Md (blue) never enters the accretion 
regime and its period converges to ~ 0.5 s. The model source with 
the highest Md (green) enters the accretion regime early on and 
after 10 5 years it has P > 30 s. The model represented by the 
solid curve can reproduce the properties of PSR J1734— 3333 (see 
Figure 3). The rectangle shows the current position of the source 
on the P — P diagram. The size of the rectangle represents the 
uncertainty in the measured P value. 



2006), the disk mass Ma, and the dipole magnetic field 
strength Bo on the surface of the star. In our earlier work 
and in the present paper, we take T p ~ 100 — 200 K. Since T p 
together with C determine the end of the active disk phase, 
they have almost no effect on the evolution during the initial 
radio phase. The initial period can affect the time of onset 
of accretion, but similar long term evolutionary curves could 
be found for a large range of initial periods. To follow the 
evolution of model sources in the initial radio phase, there 
remain only two important parameters: Md and Po- We re- 
peat the simulations tracing Md and Bo to find the entire 
allowed ranges of these parameters that can produce the 
properties of PSR J1734— 3333 simultaneously. 

Our results indicate that PSR J1734— 3333 could be at 
an age in the range ~ 3 x 10 4 — 2 x 10 5 yr depending on the 
actual initial conditions of the star. For Bo — 2 x 10 12 G and 
Po = 300 ms, M d = 3 x 10" 7 M Q gives the observed proper- 
ties of PSR J1734— 3333 simultaneously, while for Po = 100 
ms the disk mass is Md = 7.3 x 10 -7 Mq and for Po = 50 
ms Af d = 8.5 x 1O" 7 M . For B = 2 X 10 12 G and P = 300 
ms, illustrative model curves obtained with different disk 
masses are given in Figure 3. The solid (red) curve is the 
same as the tracks given in Figures 1 and 2 with the same 
Po and Po. It is seen that the model with the greatest Md 
starts the accretion phase earlier, while the source with the 
smallest Md never enters the accretion phase. These mod- 
els, represented by dashed and dot-dashed curves in Figure 
3, cannot produce all present properties of PSR J1734— 3333 
simultaneously. 

In later phases of evolution, radio pulsars similar to 



PSR J1734— 3333 are likely to be observed either as tran- 
sient AXPs in the accretion phase, if they are detected by 
means of soft gamma-ray bursts, or as XDINs after the ter- 
mination of the accretion phase, while the disk torques are 
still effective. We also present the evolution of these model 
sources on the P — P diagram in Figu re 5. A paper on the 
evolution of XDINs is in preparation (|Ertan et al.ll2013l ). 

Our long-term evolution model of neutron stars 
with fallback disks diff ers from earlier f allbac k disk 
models described in Chatteriee et alj |2000l ) and 
IMenou. Perna fc Hernquist ( 2001al ). as well as lAlparl 
|200ll ) and lAlpar et alj (|200ll ) in two crucial ingredients of 
the models: (1) the minimum temperature (T p ) at which 
the viscous activity stops and the disk becomes passive, 
and (2) the torque expression employed in the models, 
which determine the M dependence and spin evolution of 
the neutron star. 

Accordin g to th e fallback disk model proposed by 
IMenou et ail (|2001ah . the disk becomes inactive at ac- 
cretion rates at which the disk cools through the ther- 
mal viscous instability. To explain the braking indices 
of young neutron stars w ith hig h X-ray luminosities, 
IMenou. Perna fc Hernquist] (|2001bl ) employed the same 
model. In this model, the minimum mass-flow rate for the 
disk to remain active and interact with the magnetic field 
is a few times 10 14 g s _1 . The properties of the tran- 
sient AXPs, discovered later with inferred quiescent accre- 
tion rates m u ch low er than the minimum critical limit of 
IMenou et ail (|2001ah . seem hard to be explained together 
with persistent sources in a single coherent picture. An- 
other difficulty arises when applying this model to PSR 
J1734— 3333 and to other pulsars with low luminosities. At 
present, the inner disk of PSR J1734— 3333 cannot enter 
the light cylinder, and therefore the pulsed radio emission 
is allowed, provided that the mass-flow rate of the disk is 
less than about 10 13 g s _1 (even for a relatively low dipole 
field strength of 10 12 G). With this mass-flow rate, the disk 
is not active i n this model. To put it in other words, the 
model used in IMenou et alj l|2001al ) cannot be applied to 
PSR J1734— 3333. A similar problem e merge s in the fallback 
disk model used in IChatteriee et alj ([2000), since the disk 
is assumed to make a transition to ADAF regime at mass 
inflow r ates even higher than the minimum critica l rates 
given in IMenou et all (|2001al ). In IChatteriee et all (|2000h 
the disk is assumed to have no contribution to the torque 
in the radio phase, that is, this model cannot be applied to 
PSR J1734— 3333 either. In our model we do not encounter 
this problem: the disk becomes passive at low temperatures 
P p ~ 100 — 200 K, in a ccordance with theoretical work by 



llnutsuka fc Sand (|2005l), starting from th e out ermost disk 
as des cribed in detail in lErtan et alj (|2009h and lAlpar et al.l 
(20111) . 

The evolutionary model curves of the rotational prop- 
erties of neutron stars and their dependence on mass flow 
rate of the disk are sensitive to the disk torque models 
employed in the calculations. The motivation leading us 
to the particular torque model employed in earlier work 
l|Ertan fc Erkutll200Sl ; Wan et all 12003 : lAlpar et al.ll201lh 
and in the present work was the peculiarity of the con- 
temporaneous X-ray luminosity and period evolution of 
the transient AXP XTE J1810-197 in the X-ray enhance- 
ment/outburst phase. Nevertheless, the observations of XTE 
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J1810-197 l|Gotthelf fc Halpernl2007h show that in the range 
of accretion rates from a few 10 13 to 10 15 g s _1 , the torque 
remains almost independent of mass-inflow rate. We use the 
same torque model that is in good agreement with this ob- 
servation, in all subsequent work on spin evolution of neu- 
tron stars with fallback disks. This behaviour, observed in 
XTE J1810-197, cannot be accounted for by earlier disk 
torque models. 



4 CONCLUSIONS 

We have shown through detailed analysis that the observed 
period, the first and second period derivatives, and the X- 
ray luminosity of PSR J1734— 3333 can be simultaneously 
reached by a neutron star evolving with a fallback disk. The 
model is compatible with the pulsed radio emission of the 
source, since the present source properties are reached at 
a time when the accretion of matter from the disk is not 
allowed yet. 

The observed properties of PSR J1734— 3333 can be ob- 
tained at an age of ~ 3 x 10 4 - 2 x 10 5 yr, with a range of 
initial periods and with disk masses (9 x 1O _8 M0 < Md < 
1.5 x 1O _6 M0), relatively low compared to those estimated 
for AXPs/SGRs. Dipole fields with strength in the range 
9 x 10 11 — 1 x 10 13 G, together with appropriate disk masses 
and initial periods, give reasonable model curves. 

PSR J1734— 3333 has not been observed in IR to the 
best of our knowledge. A search for IR emission from the 
fallback disk in PSR J1734— 3333 will be of great interest. 
The infrared emission from the disk depends on several fac- 
tors, like inclination angle of the disk, irradiation efficiency, 
outer and inner disk radii. We will present elsewhere the 
expected IR emission from PSR J1734— 3333 and sources 
similar to it in detail, together with a discussion of the as- 
sumed parameters. 

We expect that PSR J1734— 3333 will evolve into the ac- 
cretion phase within another ~ 10 4 yr (see Figures 1 and 2). 
By the onset of the accretion phase, the mass- flow rate of the 
inner disk will have decreased to a very low level such that 
accretion does not significantly contribute to the total X- 
ray luminosity. Radio pulsars following evolutionary curves 
similar to that of PSR J1734— 3333 could be detected in 
the accretion phase if they show soft gamma-ray bursts like 
AXPs/SGRs - these would be identified as transient AXPs. 
Such sources are not likely to emit radio waves, having spun 
down to long periods by this late stage of evolution. Only 
if they are sufficiently close to Earth we could detect and 
identify them as XDINs. Some of these sources could evolve 
to pe riods longer than the AXP/SGR periods (|Ertan et al.l 
120091 ). It is hard to detect them due to low X-ray luminosi- 
ties. In short, sources like PSR J1734— 3333 may evolve to 
become transient AXPs or XDINS. The evolution of XDINs 
and their connection with AXPs/SGRs will be the subject 
of our future work. 
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